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We construct an approximate renormalization scheme 
for Hamiltonian systems with two degrees of freedom. 
This scheme is a combination of Kolmogorov-Arnold-Moser 
(KAM) theory and renormalization-group techniques. It 
makes the connection between the approximate renormaliza- 
tion procedure derived by Escande and Doveil, and a system- 
atic expansion of the transformation. In particular, we show 
that the two main approximations, consisting in keeping only 
the quadratic terms in the actions and the two main reso- 
nances, keep the essential information on the threshold of the 
breakup of invariant tori. 

PACS numbers: 05.45. -fb, 64.60. Ak 



I. INTRODUCTION 

In 1981, Escande and Dovcil |l| set up an approximate 
renormalization scheme for Hamiltonian systems with 
two degrees of freedom, in order to study the breakup of 
invariant tori, and especially to compute the threshold of 
stochasticity. Their scheme was motivated by Chirikov's 
resonance overlap criterion and by Greene's results ||] 
about the link between the existence of a torus with 
the stability of neighboring periodic orbits. They es- 
tablished the relevance of a sequence of these periodic 
orbits for the breakup of invariant tori, by setting an 
approximate transformation which focuses successively 
on smaller scales, i.e. acting like a microscope in phase 
space. 

Due to the complexity of the phase space of a non- 
integrable Hamiltonian, their method requires strong ap- 
proximations to obtain explicit expressions. Basically, 
two approximations were involved: 

(1) a quadratic approximation in the actions: their trans- 
formation produces terms that are higher than quadratic 
in the actions; in order to remain in the same class of 
Hamiltonians, they neglect these higher order terms, 

(2) a two-resonance approximation: they only keep the 
two main resonances at each iteration of the transforma- 
tion. 

The idea was to keep only the most relevant features of 
the mechanism of the breakup of a given torus. 
In this article, we construct an approximate scheme us- 
ing the same two approximations. We establish the con- 
nection between Escande's scheme |lM, ^-'^d the KAM- 
RG transformation derived in Refs. ^-0. The aim is to 
show that an exact renormalization transformation can 
be approximated by a simple transformation: It can be 



useful to derive approximate explicit expressions of uni- 
versal parameters, and to see what are the most relevant 
terms responsible for the breakup of invariant tori. The 
results we obtain support the general idea that the irrele- 
vant terms of the renormalization transformation can be 
eliminated with little loss of accuracy in the parameters 
associated to the breakup of invariant tori. 
The transformation TZ we define has two main parts: a 
KAM transformation which is a canonical change of co- 
ordinates that reduces the size of the perturbation from 
e to £^ , and a renormalization transformation which is a 
combination of a shift of the resonances and a rescaling 
of momentum and energy. 

It acts on the following class of Hamiltonians with two 
degrees of freedom, quadratic in the action variables 
A — [Ai,A2), and described by three even scalar func- 
tions of the angles = ((/5i, (f2)'- 



H{A,^)^\{\ + m{^)){n-Af 



(1.1) 



where m, g, and / are of zero average. The vector 
ujo is the frequency vector of the considered torus and 
Jl = (1,0;) is some other constant vector not parallel to 
Wq- The perturbation (m, 5, /) is of order 0(e). 
The renormalization-group approach is based on the fol- 
lowing general picture: The idea is to construct the trans- 
formation 7?. as a generalized canonical change of coor- 
dinates acting on some space of Hamiltonians such that 
the iteration of TZ converges to a fixed point. If the per- 
turbation is smaller th an c ritical, TZ should converge to 
a Hamiltonian of type (LI) with {m,g,f) = 0, which is 
integrable, and the equations of motion show that the 
torus with frequency vector is located at A = 0. All 
Hamiltonians attracted by this trivial fixed point have an 
invariant torus of that frequency (this can be considered 
as an alternative version of the KAM theorem |^ ) . If the 
perturbation is larger than critical, the system does not 
have a KAM torus of the considered frequency and the 
iteration of TZ diverges. The domain of convergence to 
the trivial fixed point and the domain of divergence are 
separated by a critical surface invariant under the action 
of TZ. The main hypothesis of the renormalization-group 
approach is that there should be another nontrivial fixed 
point (or more generally, a fixed set) on this critical sur- 
face, that is attractive for Hamiltonians on that surface. 
From its existence, one can expect to deduce universal 
properties in the mechanism of the breakup of invariant 
tori. 

Many aspects of this general picture are still at the stage 
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of conjecture, supported by some results in the pertur- 
bative regime [||, by numerical works 10 1, and by 

analogies with the related problem for area-preserving 
maps [ p]jl2| . In particular, the relation between the 
properties of the nontrivial renormalization fixed point 
and the geometric properties of the invariant torus at the 
instability threshold are not well established. The coin- 
cidence of the critical coupling of one-parameter families 
at which a torus breaks up, with the boundary of attrac- 
tion of the trivial fixed point is supported by numerical 
studies 



In Sec. ^ we describe the KAM part of the transfor- 
mation, and we make explicit the two approximations in- 
volved in this scheme: the quadratic approximation and 
the two-resonance approximation. In Sec. |lT^ , we present 
the renormalization transformation which is a combina- 
tion of the KAM part, a shift of the resonances, and 
rescalings of actions and energy. In Sec. IV, we give our 



numerical results, and in particular, we show that the ap- 
proximate scheme contains the essential features of the 
exact one. 



II. KAM TRANSFORMATION 

We perform a canonical transformation Up : {if, A) i— > 
A') defined by a generating function F{A', ip) [ p^p4| 
characterized by a scalar function X of the action and 
angle variables, of the form 



F{A\^) = A'-v + X{A\v), 



leading to 



dF 
dip 
dF 
'dA' 



A' 



dX 
dip ' 
dX 
'dA' 



(2.1) 

(2.2) 
(2.3) 



The function X is constructed such that in 7i o Up the 
perturbation terms of first order in e are equal to zero. 
Inserting Eq. (2.2) into Hamiltonian (1.1), one obtains 
the expression of the Hamiltonian in the mixed represen- 
tation of new action variables and old angle variables: 

H{A', ip) = {Vl- A'f/2 + Wo • A' 



^u}{A') 



dX 
dip 



where 



u:{A') = u;o + (O • A')0, 
h{A\ip)^]^m{ip){VL-A'f 



+ h{A\ip) + 0{e^), (2.4) 
(2.5) 

-g{ip)^-A' + J{ip). (2.6) 



The equation that determines X is thus: 
dX 

uj{A')-'^ + h{A\v)^Q. 



(2.7) 



We recall that the functions m, and / are of order 0(e) ; 
as a consequence, X is also of order 0(e). Equation (2.7) 
has the solution 



X{A',ip)^ J2 ^.(A')sin(«^.(p), 



where, if we write h{A, (p) = f^i^i-^) cos(i/ • (p) 

a;(A') • v 



(2.8) 



(2.9) 



The denominator of X^, depends on the actions: thus, 
by power expansion, it generates terms that are higher 
than quadratic in the action s. In order to remain in the 
same space of Hamiltonians (1.1), we expand the Hamil- 
tonian to the second order in the actions and neglect the 
order 0{A'^). The justification for such an approxima- 
tion is that we are interested in the torus with frequency 
vector ljq which is located at A = for the trivial fixed 
point. 



We consider Hamiltonians (1.1) with only two Fourier 
modes which are the two main resonances defined as fol- 
lows: For a frequency vector ujq, the resonances are given 
by the vectors Un — {Pn,<ln) which are the sequence of 
the best rational approximations. They are characterized 
precisely by the following property: |a;o ■ Vn\ < ■ 
for any u = {p,q) ^ i^n such that \q\ < Qn+i- For the 
frequency vector = (1/7,-1) with 7 = (1 -|- \/5)/2, 
we define the two main resonances as Ui = (1,0) and 
1/2 = (1,1)- Well known properties of continued fractions 
(see, for example, |l^) imply that the vectors i^n satisfy 

■ 1 1 



the recursion relation v 



ri+l 



Nv>„ , where N 
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The scalar functions m, g and / will have expressions of 
the form 



/(V) fi^i COs(iyi • ip) + f„,_ C0S(IV2 • ip). 



(2.10) 



The KAM part of the transformation TZ will generate 
a large set of Fourier modes from these two main reso- 
nances. The purpose of the renormalization is to relate 
these two main resonances representing the main scale, 
to the next pair of resonances (called the daughter reso- 
nances) representing the next smaller scale. For the fre- 
quency vector Wo = (l/7i~l)' the daughter resonances 
are = ( 2, 1) and 1/4 = (3,2). As the canonical trans- 
formation (2.2)-(2.3) is linear in cos(i/ • ip) and sin(i>' • y), 
t^n — t'n-i+'^n-2 and the three scalar functions (m, g, /) 
are of order 0(e), the lowest order to which the reso- 
nances uy, and 1/4 are produced, is respectively O(e^) 
and 0{£^). Thus we neg lect the order ©(e") of the KAM 
transformation. 

The next step is to expre ss t he Hamiltonian in the new 
angle variables using Eq. (2^). We notice that this equa- 
tion has to be inverted. As we need the Hamiltonian ex- 
pressed in the new coordinates to order 0{e^), we have to 
invert Eq. (2.3) up to order 0(e), since the functions to be 
expressed in the new angles are already of order O(e^). 
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Moreover, these new angle variables depend on the ac- 
tions, which we develop to order 0{A'^). The develop- 
ment can thus be summarized by expressing cos[v.(fi{ip')] 
as a function of cos{i/.(p') and sin{i/.(p'), neglecting the 
orders 0{e'^,A^). 

The final step is a translation in the A' variables such 

that the linear term is again of the form u)o ■ A' . 

The Hamiltonian expressed in the new variables becomes: 

H'{A',v')^^{A + m'{^')){n-Ar 

+ {uj,+g\v')n)-A' + f'{v'), (2.11) 

where m', g' and /' are given as functions oim, g and /. 
The constant A results of the fact that the mean value 
of to' is required to be equal to zero. We notice that the 
KAM transformation does not change rj — (1, a). 
The two-resonance approximation consists in retaining 
only the two daughter resonances and neglect all the 
other Fourier modes. Therefore the approximate KAM 
transformation lAp '^s a, map acting on a low-dimensional 
space of Fourier coeflticients: 



Up (1; , g^^ , ; m^^ , g^^ , f^^ ) 



(2.12) 



The explicit expression of this map is given in the Ap- 
pendix. 



III. RENORMALIZATION TRANSFORMATION 

We construct the approximate transformation by com- 
bining two parts: a KAM transformation (to, g, /, a) i— > 
(rn',g',f',a) as defined above, and a renormalization 
(RG) consisting of a shift of the resonances and a 
rescaling of the actions and of time {m',g',f',a) t-^ 
(to", g", /", a'). The renormalization scheme described 
in this section is for a torus of frequency vector = 
(1/7,-1) where 7 = (1 -I- ^/E)/2. It is straightforward to 
adapt it to quadratic irrationals. 

The approximate KAM-RG transformation is composed 
of four steps: 

1) a KAM transformation described in the previous sec- 
tion, which is a change of coordinates that eliminates 
terms of order 0(e), where e is the size of the perturba- 
tion; this transformation produces terms of order 0(£^), 
terms that are higher than quadratic in the actions, and 
a large set of Fourier modes. We neglect terms of or- 
der 0{A'^, £^), and also, all the Fourier modes except the 
two daughter resonances 1/3 and 1/4. We notice that this 
transformation does not change fl. 

2) a shift of the resonances: a canonical change of coordi- 
nates that maps the pair of daughter resonances {v>3, V4) 
into the two main resonances (I'l, 1/2)- 

3) a rescaling of energy (or equivalently of time). 

4) a rescaling of the action variables (which is a general- 
ized canonical transformation). 



The aim of this transformation is to treat one scale at the 
time. The steps 2), 3) and 4) are implemented as follows: 
The two main resonances (1,0) and (1,1) are replaced 
by the next pair of daughter resonances (2, 1) and (3, 2), 
i.e. we require that cos[(2, 1) • tp'] — cos[(l,0) • tp"] and 
cos[(3, 2) • (f'] = cos[(l, 1) • (fi"]. This change is done via 
a canonical transformation {A',(p') 1-^ {N~'^A',N'^(p') 
with 



2 1 
1 1 



This linear transformation multiplies loq by 7^^ (since ljq 
is an eigenvector of N); therefore we rescale the energy 
by a factor 7^ in order to keep the frequency fixed at Wq. 
A consequence of the shift of the resonances is that ft is 
changed into fl' = (1, a'), where a' = {a + 1)/(q; -I- 2). 
Then we perform a rescaling of the action variables: we 
change the Hamiltonian H' into 

A' 



H'{A',<p') = XH' [—,if 



with A such that Hamiltonian (2.11) becomes of the form 
(1.1). Since the rescaling of energy and the shift iV^ 



transform the quadratic term of the Hamiltonian into 
7^(2 + a)^[A + m'{ip')]{n' ■ A'f/2, this condition leads 
to A = 7^(2-|-a)^A. This condition has the following geo- 
metric interpretation in terms of self-similarity of the res- 
onances close to the invariant torus: the rescaling magni- 
fies the size of the daughter resonances, and places them 
approximately at the location of the original main reso- 
nances 0. 

In summary, the renormalization rescales to, (7, / and 
O = (1, a) into 



to"(v) = 



to' (iv-V) 



A 



g"{^)^-i\2 + a)g' (iV^V) , 
, l + a 



(3.1) 
(3.2) 
(3.3) 

(3.4) 



The iteration of the transformation (^^) converges to 
ck^ = 7^^. It means that $7 converges under successive 
iterations to fJ* = (1,1/7), which is orthogonal to loq 
and is the unstable eigenvector of iV^ with the largest 
eigenvalue 7^. 



IV. DETERMINATION OF THE CRITICAL 
COUPLING; AN APPROXIMATE NONTRIVIAL 
FIXED POINT 

We start with the same initial Hamiltonian as in Refs. 



H{A,^) = -{n-Af+uj^-A + ef{^), (4.1) 
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where $7 = (1,0), ljq 
a perturbation 



(l/7,-l),7= (l + \/5)/2, and 



/(tp) ^ COs(l^i • if) + COs(l^2 • <P), 



(4.2) 



where Vi = (1,0) and = (1,1)- 

We take successively larger coupling e, and determine 
whether the approximate KAM-RG iteration converges 
to a Hamiltonian with {m,g,f) = 0, or whether it di- 
verges {m,g,f) — > 00. By a bisection procedure, we de- 
termine the critical coupling Ec = 0.02885. As a compari- 
son, Escande's scheme J |l gives Sc = 0.02908. The KAM- 
RG transformation [||j7Q yields Sc — 0.02759 (Greene's 
criterion gives also Ec = 0.02759). The two-resonance 
scheme gives the result within 5%. Thus, the approxi- 
mate scheme gives a fairly good description of the critical 
surface of the breakup of the torus. 

The renormalization operator has two fixed points: a 
trivial fixed point which corresponds to the Hamiltonian 
H{A, if) ~ (il, • AY /2 -I- Wo • ^, and a nontrivial fixed 
point which lies on the boundary of the basin of at- 
traction of the trivial fixed point. The critical surface 
which is the stable manifold of the nontrivial fixed point 
is of codimension 1. The relevant critical exponent is 
5 = 2.7135. This value is close to the one obtained by 
Escande et al. ^ 5 ^ 2.7480, or by KAM-RG schemes 
P,p^ 5 — 2.6502, and by MacKay for area-preserving 
maps 111 ^ = 2.6502. 

The fact that 5 is close to 7^ can be understood by the fol- 
lowing heuristic arguments: To the first main resonance 
Vi corresponds a Fourier component Af exp[i(l,0) • <p], 
and to the second one Pexp[j(l, 1) • ¥>]. The first daugh- 
ter resonance is represented by M'exp[i(2,l) • <^]. 
The transformation is a polynomial change of coordinates 
that generates a set of Fourier modes from the two main 
resonances. The way to generate the first daughter res- 
onance at the lowest order is to combine one resonance 
Ui and one resonance 1^2. For the second daughter reso- 
nance, the change of coordinates must combine one res- 
onance Ui and two resonances V2. We rescale the phase 
space in such a way that the daughter resonances be- 
come of the same size as the main resonances. These 
arguments give the following renormalization relations: 



M' 
P' -- 



-- kiMP, 



(4.3) 
(4.4) 



where fci and fc2 are two constants that depend on how 
the transformation is performed. 

An analysis of this scheme shows that it has two fixed 
points: a trivial one M = 0, P = and a nontrivial one 



The nontrivial fixed point has a 



stable manifold of codimension 1 characterized by a rel- 
evant critical exponent; the only eigenvalue greater than 
one of the linearized map at the nontrivial fixed point 
is (5 = 7^. This relevant critical exponent does not de- 
pend on fci and ^2. This is in agreement with the general 
ideas of the renormalization group. Also, it shows that 



the relevant critical exponent of the approximate and ex- 
act KAM-RG transformation should be expected to be 
close to 7^. 

For the scaling factor at the nontrivial fixed point, we 
obtain numerically A* = 19.1248. This value can be com- 
pared with A* = 18.8282 obtained in Rcfs. 018[0,|,|l§. 
Table 1 lists some of the universal parameters associated 
with the breakup of golden tori. 

The approximate renormalization transformation has an- 
other fixed set which is a cycle of period three as it had 
also been encountered in area-preserving maps 19 
and in the KAM-RG transformation |^]. This cycle is 
simply related to the nontrivial fixed point by symme- 
tries. In particular, it belongs to the same universality 
class as the fixed point. 
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APPENDIX: TWO-RESONANCE SCHEME J J _ ^ (i , ' 



FORMULAE ' ^ 2 



^vi^y^ = o ('^^i+^'s + ^i-i/a) ' (4-11) 



We denote Ci, = cos,{v ■ ip), — cos,{v ■ ip'), and s'^ = 
sin(f • if'). The expression of the Hamiltonian in the 
mixed representation of new actions and old angles is 

H{A',^) = ]^{Sl- A'f A' 



where 



! )Ci/iCj/2Ci/3 , (4-5) 



c' c' c' 

^1/1 ^2^3 



^2^1 +2^2 +1^3 



-1/1+1/2— 1^3 ^ 
1/1+1/2+2^3 



-1/1-1/2+1/3 



-"1/1—1/2— 1/3 



-1/1-1/2+1/3 

c' 

^-1/1-1/2-1/3 



(4.12) 



(4.13) 



In the foUowing formulas, the subscripts 1, 2, 3 and 4 
denote respectively the reso nanc es (1,0), (1,1), (2,1), 
and (3,2). The Hamiltonian (4.9) becomes 



P,,,,{A) = ^ ■ J^iX,,{A)[n ■ U2X,,{A)/2 

+g^.,+m^,,Q,- A], (4.6) 

V i/i 1/2 1/3 

{A) = n- uifl ■ {A)X,, {A)m„j2. (4.7) 

We recall that X is of order 0(e), and is given by Eq. 
( ^.9[ ); as a consequence, P is of order 0{e^)^ and Q is of 
order 0{e^). 

Next, the expression of this Hami lton ian in the new an- 
gles requires the inversion of Eq. (2^) to the order 0(e). 
The expression of Cv as a function of s'^ and c'^ is 



1 



where R^^v^i-^) = 1^1 



HI 

dX^^ 



(4.8) 



OA 



The Hamiltonian (1.1) expressed in the new variables 
becomes 

H'{A!,^') = ]^{n-A'f+uja-A' 



l'l,I'2 

+ E 

+ E 

1^1,1^2.1/3 



1/1 1/2 1/3 ^1/1 ^2/2 ^2/3 



-P1/2I/1 ) ^2/2 2/3 1/2 '^1/3 

(4.9) 



i/'(A',¥'')-77("-^' 



1 



^'+0(^1+^22) 



(A2+P2l)4 



(Q221+Q212+(3l22-5'i22X, (4.14) 



where S'122 = 2P22i?2i + (f2i+fi2)(i?i2 + i?22)- Equation 
(4.14) gives the expression of the Fourier coefficients of 
the daughter resonances 



A = 1 



TO 



TO,, ~ — 



p(2) 
(2) 

12 " 



P. 



(2) 



22 1 

(2) 



21 ' 



7 1/3 



5^/4 



= ; p, 



:,(2) 

;22i 

,(1) , 



'3212 



p(l) 



f = 

f 



3(0) 

12 



I («, 



p(0) 
^21 



^(0) 
J212 



)(2) _ <;.(2) 
-122 '-'122 



^122 "^122 I ) 



)(0) _ r-(O) 
'122 "-"122 



(4.15) 
(4.16) 

(4.17) 
(4.18) 
(4.19) 
(4.20) 
(4.21) 



The next approximation is the Taylor expansion of 
H'{A', ip') to the second order in the actions e.g. 

P.,., (A) = PiJ-)^ + pW^f2 . A + Pl^l in ■ Af + 0(A3). 

(4.10) 



In the next step, we neglect all the resonances different 
from the daughter resonances (2, 1) and (3, 2) using the 
following relations 
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TABLE I. Universal parameters associated with the breakup of the golden mean torus. 



Name 



Value 



Approximate scheme 



Simple scheme 
7^ fa 2.6180 
7-2 ^ 0.3820 
7^ ^ 17.9443 
1 

J-' 



Unstable eigenvalue 
Largest stable eigenvalue 

Area multiplier A* 
Inverse mass multiplier A, 
Time multiplier 



''given in Refs. |12|,| 
''derived in Refs. 



2.6502 
0.3731 
18.8282 
1.0493 



2.7135 
0.4087 
19.1248 
1.0658 

J'' 
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